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Abstract. We demonstrate that a spherical accretion onto astrophysical black 
holes, under the influence of Newtonian or various post-Newtonian pseudo-Schwarzschild 
gravitational potentials, may constitute a concrete example of classical analogue 
gravity naturally found in the Universe. We analytically calculate the correspond- 
ing analogue Hawking temperature as a function of the minimum number of phys- 
ical parameters governing the accretion flow. We study both the polytropic and 
the isothermal accretion. We show that unlike in a general relativistic spherical 
accretion, analogue white hole solutions can never be obtained in such post- 
Newtonian systems. We also show that an isothermal spherical accretion is a 
remarkably simple example in which the only one information — the temperature 
of the fluid, is sufficient to completely describe an analogue gravity system. For 
both types of accretion, the analogue Hawking temperature may become higher 
than the usual Hawking temperature. However, the analogue Hawking tempera- 
ture for accreting astrophysical black holes is considerably lower compared with 
the temperature of the accreting fluid. 



1 Introduction 

1.1 Analogue gravity and analogue Hawking temperature 

Classical black-hole analogues are dynamical fluid analogues of general rel- 
ativistic black holes. Such analogue gravity effects may be observed when 
acoustic perturbations propagate through a classical transonic fluid. At the 
transonic point of the fluid flow the so-called acoustic horizon is formed. 
The acoustic horizon resembles a black-hole event horizon in many ways. 
In particular, the acoustic horizon emits a quasi-thermal phonon spectrum 
similar to Hawking radiation. 

In his pioneering work Unruh (1981) showed that the scalar field repre- 
senting acoustic perturbations in a transonic barotropic irrotational fluid, 
satisfied a differential equation of the same form as the Klein-Gordon equa- 
tion for the massless scalar field propagating in curved space-time with a 



2 Dasgupta et al. 



metric that closely resembles the Schwarzschild metric near the horizon. 
The acoustic propagation through a transonic fluid forms an analogue event 
horizon located at the transonic point. Acoustic waves with a quasi-thermal 
spectrum will be emitted from the acoustic horizon and the temperature 
of such acoustic emission may be calculated as (Unruh 1981) 
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where k-Q is Boltzmann's constant, a s the speed of sound, ui the component 
of the flow velocity normal to the acoustic horizon and d/dn represents the 
normal derivative. The subscript 'ah' denotes that the quantity should be 
evaluated at the acoustic horizon. The temperature T a h defined by (1) is 
an acoustic analogue of the usual Hawking temperature T^: 

T ^ (2) 

H 8Tik B GM hh V ' 

and hence T a h is referred to as the analogue Hawking temperature or simply 
analogue temperature. Similarly, the radiation from the acoustic (analogue) 
black hole is dubbed analogue Hawking radiation. Note that the sound 
speed in equation (1) in Unruh's original treatment was assumed constant 
in space. 

Unruh's work was followed by other important papers (Jacobson 1991, 
Unruh 1995, Visser 1998, Jacobson 1999, Bilic 1999). As shown by Visser 
(1998), the equation of motion for the acoustic disturbance in a barotropic, 
inviscid fluid is identical to the d'Alembertian equation of motion for a 
minimally coupled massless scalar field propagating in (3+1) Lorenzian 
geometry. In particular, the acoustic metric for a point sink was shown 
to be conformally related to the Painleve-Gullstrand-Lemaitre form of the 
Schwarzschild metric. The corresponding expression for the analogue tem- 
perature was shown to be (Visser 1998) 
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Hence, to determine the analogue Hawking temperature of a classical ana- 
logue system, one needs to know the location of the acoustic horizon and 
the values of the speed of sound, the fluid velocity and their space gradients 
on the acoustic horizon. 

In the analogue gravity systems discussed above, the fluid flows in flat 
Minkowski space, whereas the sound wave propagating through the non- 
relativistic fluid is coupled to a curved pseudo-Riemannian metric. Phonons 
are null geodesies, which generate a null surface — the acoustic horizon. In- 
troduction of viscosity may destroy Lorenz invariance and hence the acous- 
tic analogue is best studied in a vorticity free dissipationlcss fluid. 
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1.2 Motivation: Transonic black-hole accretion as an analogue 
system 

Astrophysical black holes are the final stage of gravitational collapse of 
massive celestial objects. Astrophysical black holes may be roughly classi- 
fied in two categories: stellar mass black holes having a mass of a few M© 
and supermassive black holes of mass of the order of 1O 6 M or more, where 
Mq denotes the solar mass. Both the stellar mass and the supermassive 
astrophysical black holes accrete matter from the surroundings. Depending 
on the intrinsic angular momentum of the accreting material, either the 
spherically symmetric (with zero angular momentum) or the axisymmetric 
(with non-zero angular momentum) geometry is invoked to study accreting 
black- hole systems (for details see, e.g. Frank, King & Raine 1992). 

Denoting by u(r) and a s (r) the dynamical velocity and the local sound 
speed of the accreting fluid moving along a space curve parameterized by 
r, the local Mach number of the fluid is defined as M(r) = u(r)/a s (r). The 
flow is locally subsonic or supersonic if M < 1 or M > 1, respectively. The 
flow is said to be transonic if there exists at least one point r^ such that 
M(rh) = 1. At this point a subsonic to a supersonic transition takes place 
either continuously or discontinuously. The point where this crossing takes 
place continuously is called sonic point, and if this crossing is not continu- 
ous, the transition point is called shock or discontinuity. As a consequence 
of the inner boundary conditions imposed by the event horizon, the accre- 
tion onto black holes is generally transonic. The investigation of accretion 
processes onto celestial objects was initiated by Hoyle & Littleton (1939) 
who computed the rate at which pressureless matter would be captured by 
a moving star. Subsequently, a theory of stationary, spherically symmetric 
transonic hydrodynamic accretion onto a gravitating astrophysical object 
at rest was formulated by Bondi (1952) using a purely Newtonian approach 
that includes pressure effects of the accreting material. 

Since the publication of the seminal paper by Bondi in 1952, the tran- 
sonic behaviour of the accreting fluid onto compact astrophysical objects 
has been extensively studied in the astrophysics community. Similarly the 
pioneering work by Unruh in 1981 initiated a substantial number of works 
in the theory of analogue Hawking radiation with diverse fields of appli- 
cation (for a review, see Novello, Visser & Volovik 2002). Despite the fact 
that the accreting black hole can be considered a natural example of an 
analogue gravity system, an attempt to bridge these two fields of research 
has been made only recently (Das 2004). Since both the theory of tran- 
sonic astrophysical accretion and the theory of analogue Hawking radia- 
tion are related by the same physical phenomena, it is important to study 
the analogue Hawking radiation in a transonic accretion onto astrophysical 
black holes and to compute T a h- An accreting black-hole system of classical 
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analogue gravity is a unique example of classical analogue gravity which 
exhibits both the black-hole event horizon and the analog acoustic horizon. 
Hence, an accreting astrophysical black hole may be considered an ideal 
candidate to theoretically study these two different types of horizons and 
to compare their properties. 

1.3 Pseudo-Schwarzschild black-hole accretion 

Since the relativistic effects play an important role in the regions close to 
the accreting black hole, a purely Newtonian gravitational potential in the 
form $n = —GM/r is certainly not a realistic choice to describe transonic 
black-hole accretion in general. However, a rigorous investigation of tran- 
sonic black-hole accretion (even with spherical symmetry) may be quite 
complicated in a general relativistic treatment. In order to compromise be- 
tween a relatively easy handling of the Newtonian description of gravity and 
realistic but complicated general relativistic calculations, a series of 'mod- 
ified' Newtonian potentials have been introduced. The following potentials 
have been proposed: 



where the units G = c = Afbh = 1 have been used. Unless otherwise stated, 
hereafter we use these units. The radial distance r in (4) is measured in 
units of the Schwarzschild gravitational radius r g = 2GM bh /c 2 , where M bh 
is the mass of the black hole, G Newton's gravitational constant and c is the 
speed of light. The potential <?i is proposed by Paczyhski and Wiita (1980), 
<1>2 by Nowak and Wagoner (1991), #3 and <?4 by Artemova, Bjornsson & 
Novikov (1996). 

The use of these potentials retains most of the features of space-time 
around a compact object up to a reasonably small distance from the black- 
hole event horizon, and helps in reproducing some crucial properties of 
the relativistic accretion solutions with high accuracy. Hence, these poten- 
tials might be called post-Newtonian pseudo-Schwarzschild potentials (for 
details, see Artemova, Bjornsson & Novikov 1996 and Das 2002). These 
potentials were originally introduced to study the axisymmetric accretion 
(accretion discs). However, it was shown (Das & Sarkar 2001) that some of 
these potentials could efficiently describe a spherically symmetric accretion 
as far as comparison with the full general relativistic accretion is concerned. 

In this paper we study the analogue gravity phenomena in the spher- 
ical accretion onto astrophysical black holes under the influence of vari- 
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ous post-Newtonian pseudo-Schwarzschild potentials described above. We 
use the expressions 'post-Newtonian' and 'pseudo-Schwarzschild' synony- 
mously. Our main goal is to provide a self-consistent calculation of the 
analogue horizon temperature T a h in terms of the minimum number of 
physical accretion parameters, and to study the dependence of T a h on var- 
ious flow properties. In most practical situations, a complete general rela- 
tivistic treatment turns out to be almost non-tractable. In those few simple 
cases where the general relativistic treatment is easy, one can repeat the 
calculations using various pseudo-potentials, and can compare the results 
with those obtained in the general relativistic calculations. This comparison 
may help to handle more complicated accretion systems without resorting 
to the full general relativistic treatment. 

In Section 2 we calculate the location of the acoustic horizon and eval- 
uate the relevant accretion variables on the horizon as functions of the 
fundamental parameters governing the flow. Hereafter, we use the expres- 
sions 'acoustic horizon' and 'analogue horizon' synonymously. We study 
both the polytropic (adiabatic) and the isothermal accretion. In Section 3 
we derive an analytic expression for T a h as a function of various accretion 
parameters. In Section 4 we present our results numerically and in Section 
5 we give a discussion and conclusions. 



2 Calculation of various sonic quantities 

The non-relativistic equation of motion for spherically accreting matter in 
a gravitational potential denoted by # may be written as 

du du 1 dp d<P 
dt dr p dr dr 

where u, p and p are the velocity, pressure and density of the fluid, re- 
spectively. The first term in equation (5) is the Eulerian time derivative of 
the dynamical velocity, the second term is the 'advective' term, the third 
term is the momentum deposition due to the pressure gradient and the last 
term is the gravitational force. Another equation necessary to describe the 
motion of the fluid is the continuity equation 

To integrate the above set of equations, one also needs the equation of state 
that specifies the intrinsic properties of the fluid. As mentioned earlier, 
we will study accretion described by either a polytropic or an isothermal 
equation of state. 
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2.1 Polytropic accretion 

We employ a polytropic equation of state of the form 

p = Kp\ (7) 

where the polytropic index is equal to the ratio of two specific heats, 7 = 
Cp/cy- The sound speed a s is defined by 



2 d P 



■'7 < 8 » 



where the subscript s denotes that the derivative is taken while keeping the 
entropy per particle s = S/n fixed, i.e. for an isentropic process. Assuming 
stationarity of the flow, we find the following conservation equations: 

1) Conservation of energy implies constancy of the specific energy £ 

S = ^ + ; fL + *. (9) 

2) Conservation of the baryon number implies constancy of the accretion 
rate M 

M = 4tt pur 2 . (10) 

Equation (9) is obtained from (5) with (7) and (8) and equation (10) follows 
directly from (6). 

Substituting p in terms of a s and differentiating (10) with respect to r, 
we obtain 

where ' denotes the derivative with respect to r. Next wc differentiate 
equation (9) and eliminating a' s with the help of (11) we obtain 

u = -— . (12) 

u — a% ju 

If the denominator of equation (12) vanishes at a particular radial distance 
rh, the numerator must also vanish at to maintain the continuity of the 
flow. One thus finds the sonic-point condition as 



u h = a sh = W , (13) 

where r/, is the sonic point and the sphere of radius is the acoustic 
horizon. Hereafter, the subscript h indicates that a particular quantity is 
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evaluated at r/j . The location of the acoustic horizon is obtained by solving 
the algebraic equation 

*-i(^)^-** = °- ( 14 ) 

The derivative u' h at the corresponding sonic point is obtained by solving 
the quadratic equation 

(1+7) «) 2 + 2( 7 -l) 

+ (2 7 -l)^+^=0, (15) 

' h 

which follows from (12) in the limit r^r^ evaluated with the help of 
FHospital's rule. 

Finally, the gradient of the sound speed at the acoustic horizon is ob- 
tained by substituting u' h obtained from (15) into equation (11) at the 
acoustic horizon 




Evidently, for both the Newtonian flow with = —l/r and the post- 
Newtonian flow, the location of the sonic point is identical to the location 
of the acoustic horizon due to the assumption of stationarity and spherical 
symmetry. 




2.2 Isothermal Accretion 

We employ the isothermal equation of state of the form 

RT 9 , ^ 

P= P = c 2 s P, (17) 

M 

where T is the temperature, R and n are the universal gas constant and 
the mean molecular weight, respectively. The quantity c s is the isothermal 
sound speed defined by 



2 _ 

Cs dp 



= 0T, (18) 

T 



where the derivative is taken at fixed temperature and the constant — 
Kb/ipmn) with m# ~ m p being the mass of the hydrogen atom. In our 
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model we assume that the accreting matter is predominantly hydrogen, 
hence fi ~ 1. Now, the specific energy equation takes the form 

u 2 

£ = — + OT\np + $, (19) 

whereas the accretion rate is given by (10) as before. 

The radial change rate of the dynamical velocity is again given by (12). 
From equation (12) and with (18) we find the sonic point condition as 




The derivative of u is obtained from (12) by making use of l'Hospital's rule 
as before. We find 

< = -J- Ifa+n) , (2D 

where the minus sign in front of the square root indicates accretion (the 
plus would correspond to a wind solution). Note that the quantities in 
equations (20) and (21) are functions of the fluid temperature T only. 
Hence the isothermal spherical accretion can be essentially described as a 
one-parameter solution of the hydrodynamical equations, parameterized by 
T. 



3 Calculation of Analogue Temperature 

To calculate T a h for various <£, we first write down the general expression 
for T a h in the relativistic form and then we reduce the expression by taking 
the non-rclativistic post-Newtonian limit. 

The relativistic acoustic metric tensor G M „ is defined as (Moncrief 1980, 
Bilic 1999) 

= [9^ - (l - a 2 ) u^u,,] , (22) 

where n is the particle number density and h = (p + p) /n is the specific 
enthalpy. The corresponding surface gravity may be calculated as (Bilic 
1999) 



1 — at on 



(23) 

ah 



where is the stationary Killing field and d/dn is the derivative normal 
to the acoustic horizon. This expression together with the Newtonian limit 

\e\=9oo^ (l + ^j) ( 24 ) 
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gives a general expression for the temperature of the analogue Hawking 
radiation in a spherically accreting fluid in the Newtonian as well as in any 
pseudo-Schwarzschild gravitational potential 

J_ l 2c 2 + $ h 
ah ~27r Kfc V 2c 2 

The quantities required to calculate the analogue temperature (25) are 
obtained using the formalism presented in Section 2.1. The Newtonian 
and post-Newtonian polytropic spherical accretion onto astrophysical black 
holes is an example of an analogue-gravity system in which the value of T a h 
can be calculated using only two measurable physical parameters, namely, 
£ and 7. 

For a particular value of {£, 7}, we define a dimensionless quantity r as 
the ratio of T a h to Th which turns out to be independent of the black-hole 
mass. In this way we can compare the properties of the acoustic versus event 
horizon of a spherically accreting black hole of any mass. Using equations 
(11)-(16) we find 

= = 4 ! 2+<Ph ( 2 ^ fi±i] 

T T H V 2 \2-r h $ h ) \ 2 J 

A(7)-(l+ 7 )^ (26) 

y >h 

where f (7) = (0.000757 2 - 5.00157 + 3.00075). The quantities $ h , & h , and 
are obtained from equation (4) for various potentials, and 77, is calcu- 
lated from (14) for an astrophysically relevant choice of {£ ,7}. 

The quantity £ is scaled in terms of the rest mass energy. The values 
£ < correspond to negative energy accretion states where a radiative 
extraction of the rest mass energy from the fluid is required. To make such 
an extraction possible, the accreting fluid has to possess viscosity or other 
dissipative mechanisms, which may violate Lorenz invariance. The values 
£ > 1 arc mathematically allowed. However, large positive £ represent 
the flows that start from infinity with a very high thermal energy, and the 
£ > 1 accretion represents enormously hot flow configurations at very large 
distances from the black hole. These configurations are not conceivable in 
realistic astrophysical situations. Hence we set £ £ £ 1. 

Concerning the polytropic index, the value 7 = 1 corresponds to the 
isothermal accretion in which the accreting fluid remains optically thin. 
This value is the physical lower limit since 7 < 1 is not realistic in accretion 
astrophysics. On the other hand, 7 > 2 is possible only for a superdcnsc 
matter with a substantially large magnetic field and a direction-dependent 
anisotropic pressure. This requires the accreting material to be governed by 



dr 



(o s - u) 



(25) 



ah 
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magneto-hydrodynamic equations, dealing with which is beyond the scope 
of this paper. We thus set 1 ^ 7 ^ 2 for a polytropic flow, and separately 
consider the case 7 = 1, corresponding to an isothermal accretion. The 
preferred values of 7 for a realistic polytropic black hole accretion range 
from 4/3 to 5/3 (Frank et. al. 1992). 

Note that if (a' s — u')h is negative, one obtains an acoustic white-hole 
solution. Hence the condition for the existence of the acoustic white hole is 



y5 f(7) " (1+7) ^ <0 - (27) 

Since 7 and can never be negative, and since ^ and 4% are always 
real for the preferred domain of {£,7}, acoustic white- hole solutions are 
excluded in the astrophysical accretion governed by the Newtonian or post- 
Newtonian potentials. 

For an isothermal flow, the quantity c' s is zero and using (21) we find 

where r/j should be evaluated using (20). Clearly, the fluid temperature T 
completely determines the analogue Hawking temperature. Hence, a spher- 
ical isothermally accreting astrophysical black hole provides a simple sys- 
tem where analogue gravity can be theoretically studied using only one free 
parameter. 

4 Results 

4.1 Polytropic accretion 

Given a potential <£, we calculate as a function of {£, 7} using equation 
(15). The physically motivated choice of the domain {£,7} is discussed in 
Section 3. Note that the physically acceptable values of Th must be larger 
than 1 since the black-hole event horizon is located at r = 1. However, 
we show the results only for r^>2 since the pseudo-potentials (4) may not 
accurately mimic the Schwarzschild space-time too close to the black-hole 
event horizon. We find that r h anti-correlates with {£,7}, which is obvious 
from equation (14). Physically, this effect is due to a subtle competition 
between the magnitudes and radial gradients of the velocity field and the 
speed of sound. On the one hand, increase of £ at fixed 7 means increase of 
the asymptotic energy configuration of the flow which is basically thermal 
in nature. As a result, the larger the energy £, the closer to the accretor 
must the flow approach in order to gain sufficiently large dynamical velocity 
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Fig. 1. The ratio of the analogue to the Hawking temperature r = T^h/Tn as a 
function of the specific flow energy £ and the polytropic index 7 for the Newtonian 
(marked with $n) and four pseudo-Schwarzschild potentials (marked with #1, $2, 
$3 and $4). 

and its gradient. On the other hand, keeping £ fixed, increase of 7 means 
increase of the heat capacity of the accreting fluid and in turn its speed of 
sound. As a result, the flow becomes supersonic at a smaller distance to 
the black hole. 

Once we have calculated rn, we calculate the ratio r = T a h/Tn as a 
function of two flow parameters £ and 7 using (26) and plot the results in 
figure 1 for various potentials. Note that for all potentials except ^n, the 
ratio r may be larger than unity for some {£ ,7} and increases both with 
increasing £ and with increasing 7. In the region 1 < < 2, the ratio r 
may become very large, and a small zone may appear where r is again less 
than one. In the Newtonian case, r becomes larger than unity only in the 
region 1 < < 2. However, as discussed earlier, we do not rely on the 
results obtained in this region. 

In figure 2 we classify the {£, 7} parameter space to show some distinct 
important regions. The region marked with r < 1 is a range of values 
of {£,"/} for which the analogue temperature is always lower than the 
Hawking temperature. Similarly, r > 1 denotes a range of values of {£,7} 
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Fig. 2. Classification of the {£, 7} parameter space for a polytropic spherical ac- 
cretion flow in the Newtonian (marked with ^n) and in four pseudo-Schwarzschild 
potentials (marked with $1, $2, ^3 and $4). The parameter space is divided 
in four regions in which the analogue temperature Tah is always lower than 
the Hawking temperature Th (marked with r < 1), T a h is higher than Th 
(marked with r > 1), the use of the pseudo-potentials is unreliable (marked 
with 1 < r h < 2), the accretion flow is always subsonic (marked with NS). 



for which T a h is higher than Tjj. In these regions we have Th>2. The region 
marked with 1 < < 2 represents the values of {£, 7} for which a physical 
solution may be obtained for 1 < 77, < 2, and r may be larger than (mostly) 
and less than (rarely) unity. The regions marked with NS represent the 
{£,7} subset for which no physical solution r/i > 1 exists. Hence, in this 
region of the parameter space there is no transonic flow and the acoustic 
horizon does not form. Note that such a situation may occur only with ^3 
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and with purely Newtonian potentials. With all other potentials, physical 
values (r/j > 1) exist in the entire domain of {£,7} • 

The domination of the analogue temperature over the Hawking temper- 
ature is most prominent with the Paczyhski and Wiita (1980) potential, 
and least prominent with the pure Newtonian potential. If we denote by 
t* the analogue to the Hawking temperature ratio for a particular <P at a 
particular value of {£, 7} for which a physical solution rh exists for all five 
potentials, then 

r *l > T S>2 > r *3 > r *4 > r *N_ (29) 




Fig. 3. The ratio of the analogue to the Hawking temperature r = Tah/Tn as a 
function of the acoustic velocity c s (scaled in the unit of the velocity of light) at 
r h , for spherical isothermal accretion in the Newtonian (marked with #n) and 
four pseudo-Schwarzschild potentials (marked with $1, $2, &3 and $4). Results 
are shown for rn>2 only. 



4.2 Isothermal accretion 

As demonstrated in Section 2.2, the flow temperature T fully determines the 
location of the acoustic horizon rh and other relevant quantities evaluated 
at rh- We take a representative fluid temperature (in degrees of K) in the 
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range 10 8 <T<4x 10 12 , where the upper limit ensures that should never 
be less than 2. We find that the horizon radius calculated using equation 
(20) increases with T. This is obvious because the higher the temperature, 
the larger the speed of sound, and the flow needs to get closer to the 
accretor to gain the dynamical velocity u and become supersonic. For a 
given temperature T, we find 

(30) 

In figure 3 we plot the ratio r as a function of the isothermal acoustic 
velocity c s (scaled in the unit of the velocity of light) at the acoustic horizon 
r/j, for five potentials as in figure 1. The value of c s is obtained from eq. 
(20). The lines in the plot extend as far as a physical value > 2 exists. 
We find that r always increases with T (hence with c s ) and at a particular 
T (as well as at that corresponding c s ) we also find 

5 Discussion 

In this work we have demonstrated that a spherical accretion onto astro- 
physical black holes, under the influence of Newtonian and various post- 
Newtonian pseudo-Schwarzschild potentials, can be considered as an ana- 
logue gravity system naturally found in the Universe. Spherically accreting 
black holes are unique analogue gravity systems in which two kinds of 
horizons, black-hole event and acoustic analogue horizon, may be simulta- 
neously present. We have calculated the analogue horizon temperature T a h 
for spherically accreting systems as a function of the minimum number of 
physical parameters governing the accretion flow. We have modelled the 
flow in a self-consistent way and studied the dependence of T a h on various 
flow properties. In our model the calculations have been performed with a 
position- dependent speed of sound using the relativistic formalism in curved 
space (Bilic 1999). 

Note, however, that the analogy has been applied to describe the classi- 
cal perturbation of the fluid in terms of a field satisfying the wave equation 
in an effective geometry. It is not our aim to provide a formulation by which 
the phonon field generated in this system could be quantized. To accom- 
plish this task, one would need to show that the effective action for the 
acoustic perturbation is equivalent to a field theoretical action in curved 
space, and the corresponding commutation and dispersion relations should 
directly follow (see, e.g. Unruh & Schutzhold 2003). Such considerations 
would be beyond the scope of this paper. Also note that for all types of 
accretions discussed here, the analogue temperature T a h is many orders of 
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magnitude lower compared with the fluid temperature of accreting matter. 
Hence, the black-hole accretion system may not be a good candidate for 
detection of the analogue Hawking radiation. 

We have studied both the polytropic and the isothermal flow. For a 
polytropic flow, two free parameters £ and 7 arc required to calculate T a h- 
The polytropic index 7 is assumed to be constant throughout the fluid. 
A realistic model of the flow with no assumptions would perhaps require 
a variable polytropic index having a functional dependence on the radial 
distance, i.e. 7 = j(r) instead of a constant 7. However, we have performed 
the calculations for a sufficiently large range of 7 and we believe that all 
possible astrophysically relevant polytropic indices are covered. 

Das & Sarkar (2001) have shown that among all pseudo-potentials, $1 
and <?4 are the best because they provide an accurate approximation to 
the full general relativistic solution up to a reasonably close distance to the 
event horizon. While <Pi is the best approximation for an ultra-relativistic 
flow with 7 = 4/3, the potential $4 turns out to be the best approximation 
when the flow becomes non-relativistic, i.e. when 7 ~ 5/3. Hence, the 
temperature T a h calculated with <Pi and ^4 is comparable with the value of 
T a h calculated for a general relativistic flow. In the isothermal case, there 
is presently no such comparative study and hence no preference among 
various <P may be set for an isothermal accretion. 

From our calculations it follows that the main difference between the 
post-Newtonian polytropic flow (irrespective of the potential used) and 
the full relativistic flow is that an analogue white hole can never form in 
the Newtonian or post-Newtonian spherical accretion, whereas in a general 
relativistic accretion, analogue white-hole solutions exist for an extended 
parametric region of {£,7} (Das 2004). In an isothermal accretion with a 
constant speed of sound, the question of the analogue white-hole formation 
does not arise at all. 

We have found that the isothermal spherical accretion onto an astro- 
physical black hole represents a simple analogue gravity model in which 
one physical parameter, the fluid temperature T, completely specifies the 
system and provides sufficient information for the calculation of the ana- 
logue temperature. We would, however, like to raise one point of caution. 
In a realistic astrophysical situation with perturbations of various types, 
isothermality of the flow is difficult to maintain if the acoustic horizon is 
formed at a very large distance from the event horizon. Owing to the fact 
that, close to a spherically accreting Schwarzschild black hole, the electron 
number density n e falls off as r~ 3 / 2 while the photon number density n 7 
falls off as r -2 (Frank et al. 1992), the ratio of n e to n 7 is proportional to 
y/r. Hence, the number of electrons per photon decreases with decreasing 
r, so that close to the accretor, the momentum transfer by photons on the 
accreting fluid may efficiently keep the fluid temperature roughly constant. 
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Hence, isothermality may be a justified assumption at small values of rv 
However, the momentum transfer efficiency falls off with increasing r and 
the isothermality assumption may break down far away from the black 
hole. 

We have found that in both the adiabatic and the isothermal flow, the 
analogue temperature may take over the Hawking temperature, i.e. the 
ratio t may be larger than unity. This effect has also been established in 
a fully relativistic flow (Das 2004). One should note that since r decreases 
with increasing r^, high values of t are obtained if the acoustic horizon is 
formed close to the event horizon, i.e. for 1 < < 2. However, our post- 
Newtonian model may not be reliable in this region since the difference 
between a post-Newtonian metric and the general relativistic one is likely 
to be quite substantial in the vicinity of the event horizon. The potentials 
discussed here could only be used to obtain a more accurate description 
compared with the purely Newtonian approach. Hence very large values of r 
are consistent only in a complete general relativistic treatment of accretion. 

One interesting important problem is the study of the analogue gravity 
and the calculation of the analogue temperature for a non-spherical ax- 
isymmetric accretion (accretion discs). In this case, the situation is more 
complicated because, first, a realistic disc accretion requires a more in- 
volved parametric space and, second, more than one sonic points exist for 
specific boundary conditions (Das 2002, Das, Pendharkar & Mitra 2003, 
Das 2004a, Barai, Das & Wiita 2004). A detailed study of these issues is 
beyond the scope of this article and is presented elsewhere (Abraham, Bilic 
& Das). 
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